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Observational Goals for Max ’91 to Identify 
the Causative Agent for Impulsive Bursts 




D. A. Batchelor (NASA/GSFC) 

Recent studies of impulsive hard X-ray and microwave bursts suggest that a propagat- 
ing causative agent with a characteristic velocity of order 1000 km s is responsible 
for these bursts. In this presentation, the results of those studies will be summarized 
and observable distinguishing characteristics of the various possible agents will be high- 
lighted, with emphasis on key observational goals for Max ’91 campaigns. The most 
likely causative agents suggested by the evidence are shocks, thermal conduction fronts, 
and propagating modes of magnetic reconnection in flare plasmas (although other pos- 
sible agents cannot as yet be ruled out). With the new instrumentation planned for 
Max ’91, high spatial- resolution observations of hard X-ray sources have the potential 
to identify the agent by revealing detailed features of source spatial evolution. Coor- 
dinated observations with the Very Large Array and other radio imaging instruments 
are also obviously of great importance, as well as detailed modeling of coronal loop 
structures to place limits on density and temperature profiles in the loops. With the 
combined hard X-ray and microwave imaging observations, aided by loop model results, 
the simplest causative agent to rule out would be the propagating modes of magnetic re- 
connection. To fit the observational evidence, reconnection modes would need to travel 
at approximately the same velocity (the Alfven velocity) in different coronal structures 
that vary in length by a factor of 10 3 . Over such a vast range in loop lengths, it is diffi- 
cult to believe that the Alfven velocity is constant. Thermal conduction fronts would be 
suggested by sources that expand along the direction of B and exhibit relatively little 
particle precipitation. Particle acceleration due to shocks could produce more diverse 
radially expanding source geometries with precipitation at loop footpoints. 
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I. INTRODUCTION 


The main objective of research on solar and stellar flaxes is the discovery of 
the process responsible for the unpredictable, rapid, enormous releases of energy that 
occur in flares. The clearest observational clues are the impulsive X-ray and microwave 
radiations emitted during flares because these radiations offer the most direct informa- 
tion about the energy release process that boosts the radiating particles to an energy 
per particle E > 25 keV (see Fig. 1). 

There is ample evidence that simultaneous impulsive bursts of hard X rays 
(photon energy range from 25 keV to 500 keV) and microwaves (3 GHz < / < 100 GHz 
band) both radiate from one distribution of energetic electrons (Takakura and Kai 1966; 
Matzler 1978; Gary 1985; Schmahl, Kundu, and Dennis 1985; Kai 1986). The X rays 
are bremsstrahlung and the microwaves are gyrosynchrotron radiation. Crannell et 
al. (1978) made the assumption of a common source electron distribution, interpreted 
the hard X-ray spectra as thermal bremsstrahlung, and introduced a method to derive 
a characteristic length scale of a burst source. Their analysis showed that the rise 
times of solar impulsive bursts were correlated with spatial length scales of the sources. 
The correlation was confirmed by Batchelor et al. (1985) and Batchelor (1987), using 
SMM observations of hard X rays and microwave observations from Bern, Toyokawa, 
and Itapetinga observatories. The correlation is linear, consistent with the thermal 
conduction front model (TCF, see Brown, Melrose, and Spicer 1979). 

A similar method for deriving source lengths was used in the work described 
herein, but it was adapted for use with nonthermal source models. The observations 
analyzed by Batchelor et al. (1985) and Batchelor (1987) were re-analyzed, employing 
the standard nonthermal models in the literature, the thick-target model (TT, Brown 
1971) and trap- plus-precipitation model (TP, Melrose and Brown 1976). 

II. METHOD OF ANALYSIS 

The analysis method is briefly described in this section; a more detailed descrip- 
tion of the method was given by Batchelor et al. (1985) and Batchelor (1989). The hard 
X-ray observations were made with the Hard X-ray Burst Spectrometer (HXRBS, see 
Orwig, Frost, and Dennis 1980). The microwave observations were made at Bern (Ma- 
gun et al. 1981), Toyokawa (Torii et al. 19(9), and Itapetinga (Kaufmann et al. 1982). 

If the microwave spectrum of a burst is known throughout a sufficiently broad 
band, it generally exhibits a low-frequency segment with positive spectral index a (i.e. 
flux density 5 oc / Q , where / is frequency in Hz) and a high-frequency segment with 
negative a (the full spectrum is termed “C-type”, Guidice and Castelli 1975; see Fig. 2 
for examples). The positive-a segment is usually attributable to emission that is op- 
tically thick, due to self- absorption, and the negative-c* segment to emission that is 
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Figure 1. An impulsive solar burst of hard X rays and mi- 
crowaves (HXRBS and Bern data respectively). The descend- 
ing arrow marks the end of the first steady, impulsive rise 
during this burst, and the flux at this time was used to de- 
termine t r (see text). 



MICROWAVE FREQUENCY (GHz) 

Figure 2. Examples of two microwave burst spectra, each at 
the time of hard X-ray burst maximum (Bern data), (a) 1981 
August 10, (b) 1981 May 4. Arrows indicate points used to 
determine S 2 and f 2 (see text). 


134 



optically thin (Matzler 1978; Ramaty 1969). The optically-thick emission was used here 
to derive an area of each source at the time when the counting rate summed over the 
whole HXRBS energy range was maximum. 

First, 30 bursts observed with HXRBS and the radiotelescopes were selected. 
Bursts were selected with steady rises to a peak or plateau in their time histories (t r < 
30 s), and microwave spectra observed in the positive-a segment. (Twice the time from 
half-maximum to maximum was defined a s t r .) For each burst, values denoted 52 and 
/2 were selected at the highest-frequency observation on the positive-a segment (see 
Fig. 2). The X-ray spectrum from HXRBS of each burst at its peak flux was fitted with 
the power law 1(e) = A\ e~' 1 photons cm -2 s -1 keV -1 , in order to obtain the hard 
X-ray spectral index 7. 

The approximate formulae derived by Dulk and Marsh (1982) for optically- 
thick microwave flux from a power law distribution of electrons ( N(E ) oc J E~° elec- 
trons cm -3 ) were used to compute the effective emission temperature of the nonthermal 
microwave source T e ff (deg K) and A, the area of each source (cm 2 ): 


T eff = 2.2 x 10 9 10 -031 *■ 

( 1 ) 
( 2 ) 


(sin 8) 


- 0 . 36 - 0.06 6 


(///* j 


0 . 50 + 0.085 6 


A = 7.5 X 10 43 5 / 


r-2 rp— 1 

i e // 5 


In Eq. (1), 6 is the power-law index of the electron distribution. In the TT model, 
S = 7 + 1, and in the TP model, 6 = 7 — B is the magnetic field, the variable fs 
= 2.8 x 10 6 B is the gyrofrequency (Hz), and 8 is the angle between the B vector and 
the line of sight. In Eq. (2), the Sun- Earth distance is accounted for, and the units of 
5 are solar flux units (1 SFU = 10 -22 W m -2 Hz -1 ). 

The uncertainties in B and 8 were treated as follows. Typical values of B in 
flares, deduced from other solar observations, generally range from 100 to 1000 gauss 
(Svestka 1976; Brown, Smith, and Spicer 1981). In this work, two values of L = A 
were computed for each burst, with B set equal to each of those two limits; the range in 
results was treated as the ±lcr uncertainty. L is admittedly a crude approximation of 
the characteristic source dimensions, but its large range of variation - more than three 
orders of magnitude - makes it useful despite uncertainties in length-to-width ratio of 
the sources that it incorporates. The angle 8 could not be measured. Gyrosynchrotron 
radiation is emitted most intensely in directions « 90° from B (see Eq. (1)), so that 
the part of each source with 8 nearest to 90° contributes most of the radiation. The 
dependence of L on 8 is slight in the range 45° < 8 < 90°, so 8 — 45° was assumed for 
each burst, as an estimate of the mean value. Little error is incurred since bursts with 
small 8 would have low T e //, and tend not to be observed. 
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III. RESULTS 


The observational parameters and results appear in Table 1. Figure 3 shows 
the plot of L vs. t r for the TT model. The TP results make a very similar plot in 
Figure 4, except for a reduction of the L values by about 30%. Three of the flares 
occurred on or beyond the solar limb, and are marked with square symbols. Because 
these bursts might have been partially occulted by the limb, altering A , only the other 
27 were included in the least-squares fit (Bevington 1969), drawn as the straight line in 
the plot. The slope of the best fit line, 1.07 ± 0.05, is unity within the uncertainties. 

At least one of the limb bursts, the 1981 Dec. 7 event, has a value of L that 
is markedly displaced from the correlation, appearing in Fig. 3 with a square symbol 
well below the best fit line. The displacement of its value of L is in the direction 
consistent with partial occultation of a flare located just beyond the solar limb. This 
is an important test of consistency in the interpretation of A as representative of the 
true source area; displacement of a limb flare value of L in the other direction relative 
to unocculted burst values would have been an unresolvable conflict with geometry. 

As described by Batchelor et al. (1985), the parameters entering the corre- 
lation were checked to reveal possible more fundamental underlying correlations. Of 
particular concern was the “Big Flare Syndrome” (Kahler 1982), the possible associ- 
ation of larger flare intensities with larger source sizes, harder X-ray spectra, etc. No 
underlying parameter correlations were capable of accounting for the close least-squares 
fit of the relation in Fig. 3. The “Big Flare Syndrome” in particular was ruled out by 
broad scatter in the plot of S 2 vs. t r . 

IV. DISCUSSION 

The deduction that L oc t r if one chooses the TT or TP model has important 
implications for the interpretation of impulsive flare phenomena. The correlation of 
L with t r was already known on the basis of thermal flare models. The differences 
between thermal and nonthermal models for impulsive microwave and X-ray emission 
are substantial, and the (L,t r ) correlation must result in disparate ways from each of 
these models. The key inference from these results is: whether one chooses a model 
from either the nonthermal or thermal alternatives, there is clear evidence in each case 
that all impulsive bursts in the 0.1 < t r < 30 s range are due to one type of causative 
agent, characterized by a velocity of order 10 8 cm s -1 . A precise single velocity for the 
agent is not implied by the data, of course, given the crudeness of L and the scatter of 
the values; nevertheless, a type of agent that travels at speeds within less than an order 
of magnitude range is suggested by Fig. 3, and this will be termed “a characteristic 
velocity” hereafter. 
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TABLE I 


Observational and Derived Parameters 
of Solar Hard X-ray/Microwave Bursts 


Date UT y t r 

(s) 


80 

Mar 

29 

0918 

09 

3.3 

3.0 

80 

Mar 

29 

0955 

06 

2.8 

5.2 

80 

Jun 

04 

0654 

19 

3.8 

7.0 

80 

Jun 

29 

1041 

35* 

3.6 

3.6 

80 

Jul 

01 

1626 

53 

2.8 

0.8 

80 

Jul 

01 

1626 

56 

2.9 

0.9 

80 

Jul 

01 

1626 

59 

3.2 

0.8 

80 

Jul 

01 

1627 

02 

2.8 

1.0 

80 

Jul 

01 

1627 

04 

3.0 

0.8 

80 

Jul 

01 

1627 

08 

2.6 

1.0 

80 

Jul 

01 

1627 

13 

2.6 

1.4 

80 

Oct 

09 

1123 

58 

4.2 

5.2 

80 

Nov 

05 

2233 

02 

3.8 

24. 

80 

Nov 

06 

0650 

51 

4.5 

20. 

80 

Nov 

08 

1450 

25 

5.9 

7. 

80 

Nov 

18 

0718 

08* 

3.4 

2.2 

80 

Dec 

17 

0845 

37 

3.2 

3.2 

81 

Mar 

23 

0655 

49 

3.5 

6.0 

81 

Apr 

10 

1644 

53 

4.3 

10. 

81 

Apr 

15 

0643 

09 

4.7 

3.8 

81 

Apr 

18 

1049 

28 

4.5 

5.0 

81 

Apr 

26 

1115 

31 

4. 

11. 

81 

May 

04 

0838 

03 

4.2 

1.8 

81 

Jul 

19 

0533 

25 

3.9 

12 . 

81 

Jul 

20 

1311 

27 

4.5 

22. 

81 

Aug 

10 

0658 

50 

3.9 

2.6 

81 

Dec 

07 

1451 

02* 

3 . 1 

10.0 

84 

May 

21 

1326 

29 

2.7 

0 . 2 (m) 

84 

May 

21 

1326 

30 

2.4 

0 . 1 (m) 

84 

May 

21 

1326 

37 

3.2 

0 . 1 (m) 


S 2 f 2 L 

TT model TP model 


(SFU) 

510. 

(GHZ) 

10.4 

(cm) 
6.5 X 

10 8 

10 8 

4.4 

(cm) 
x 10 8 

210. 

10.4 

3.6 

X 

2.5 

X 

10 H 

320. 

8.4 

8.0 

X 

10 8 

10° 

10® 

10 7 

5.3 

X 

10 8 

31. 

8.4 

2.4 

X 

1.6 

X 

10 8 

135. 

19.6 

1.2 

X 

8.4 

X 

10 7 

460. 

19.6 

2.3 

X 

1.6 

X 

10 8 

39. 

19.6 

7.1 

X 

10 

5.0 

X 

10 7 

82. 

35.0 

4.1 

X 

;?e 

10 8 

1 0 Q 

3.0 

X 

10 7 

486. 

28.0 

1.4 

X 

1.0 

X 

10 8 

894. 

35.0 

1.3 

X 

9.5 

X 

10 7 

1330. 

35.0 

1.6 

X 

1.2 

X 

10 8 

100. 

5.2 

1.0 

X 

' 10 9 

6.5 

X 

10 8 

2400. 

9.4 

1.9 

X 

10 9 

1.2 

X 

10 9 

910. 

8.4 

1.6 

X 

10 8 

1.1 

X 

10 9 

33. 

8.4 

4.6 

X 

10 7 

3.0 

X 

10 8 

44. 

19.6 

7.9 

X 

10 8 

10® 

10 9 

5.5 

X 

10 7 

280. 

8.4 

6.3 

X 

4.2 

X 

10 8 

260. 

8.4 

6.6 

X 

4.4 

X 

10 8 

120. 

5.2 

1.1 

X 

10 8 

10 s 

7.3 

X 

10 8 

14 . 

5.2 

4.4 

X 

2.8 

X 

10 8 

55. 

8.4 

4.0 

X 

10 9 

2.7 

X 

10 s 

75. 

3.2 

1.7 

X 

10 9 

10 8 

1.0 

X 

10 9 

17. 

5.2 

4.2 

X 

2.7 

X 

10 8 

1400. 

19.6 

5.0 

X 

10 9 

3 . 5 

X 

10 s 

170. 

2.8 

3.6 

X 

10 8 

10 8 

10 f 

2.2 

X 

10 9 

200. 

8.4 

6.5 

X 

4.3 

X 

10 8 

240. 

19.6 

1.7 

X 

1.2 

X 

10 8 

20. 

90.0 

5.3 

X 

10? 

4 . 1 

X 

10 6 

50. 

90.0 

7.9 

X 

10* 

6.1 

X 

10 6 

30. 

90.0 

7.0 

X 

io 6 

5.4 

X 

10 G 


* Limb flares, excluded from length-rise time correlation due to 
possible limb occultation. 

m The 1984 May 21 bursts had rise times shorter than HXRBS time 
resolution. The microwave rise time has therefore been used. In 
all other cases, the rise time is derived from HXRBS data. 
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MICROWAVE SOURCE LENGTH SCALE (cm) 


LENGTHS VS. RISE TIMES: THICK TARGET 


i o 10 

10 9 

10 8 
10 7 

10 6 

0.01 0.1 1 10 100 

HARD X-RAY RISE TIME (s) 

Figure 3. Plot of burst source length scales L vs. rise times 
t r , assuming the TT model. Symbols: crosses label three 
bursts from the 1984 May 21 flare observed at Itapetinga; 
diamonds label seven bursts from the 1980 July 1 flare; all 
other symbols label bursts from other flares. Squares indicate 
limb events; because these bursts might have been partially 
occulted (reducing L ), they were excluded from the least- 
squares fit. The fit, performed with the function log L = 
log a + b log t r , yielded a = 10 8 cm s -1 , 6 = 1.07, oi og a = 
0.08, <T{, = 0.05, correlation coefficient r = 0.95. There are N 
= 27 points, so the probability that the quantities L and t r 
are actually uncorrelated is P c (r,N) <C 10 -6 . 
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Figure 4. Same as Fig. 3 for the trap-plus-precipitation (TP) 
model. The fit parameters are similar: a = 7.3 x 10 7 cm s 
b = 1.03, ai 0 g a = 0.07, o b = 0.04, correlation coefficient r 
— 0.95. Again with N = 27 points, the probability that the 
quantities L and t r are actually uncorrelated is P c (r,N) <C 
10 " 6 . 
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The nature of the causative agent and its characteristic velocity are specific to 
each model. The mean values are 10 8 cm s -1 for the TT model, 7 x 10 7 cm s -1 for 
the TP model, and 2 x 10 8 cm s -1 for the TCF model. In each model, the role of this 
agent is fundamental: it couples the energy release process to the geometrical scale of 
the burst source region and determines the duration of the most rapid impulsive energy 
release phase. Possible agents are shocks in the nonthermal models and TCFs. Any 
other processes of flare energy release that may be proposed also must include such an 
agent. 

Perhaps it is no great surprise that a velocity of order 10 8 cm s -1 characterizes 
the growth of flare sources, but this has not been demonstrated heretofore by means of 
standard nonthermal models, nor has the dependence of the characteristic velocity on 
the chosen model been derived before. One can envision many ways that flare sources 
could grow with some characteristic speed, so the results described above do not permit 
us to discriminate between models. Nevertheless, the evidence for interpretation of 
the burst rise time as a phase of growing source area is an important departure from 
the long-standing interpretation that the time behavior of the injection of accelerated 
particles, convolved with particle propagation effects, determines burst rise times (e.g., 
Emslie 1983; Lu and Petrosian 1988). 

The results of this work suggest other new lines of investigation. First, the 
correlation L oc t r supports the future use of A and L measurements as meaningful 
tools of source analysis, whereas they were more questionable before. Future studies 
may benefit from tests for correlations between L or A and other available flare param- 
eters. Second, the causal agent responsible for the correlation must be thought of as a 
physical process that can operate in a similar manner on size scales that vary by three 
orders of magnitude, a useful datum for investigations of simulated flare energy release 
via magnetic reconnection. Conversely, something must serve to terminate the process 
on equally wide-ranging size scales. Third, these results give information about a few 
sources that are less than one arc second in angular size. Even the most advanced mi- 
crowave imaging system for flares, the Very Large Array, can’t provide such information 
with its images. 

The clear evidence herein for source growth with a characteristic velocity dur- 
ing burst rises should be tested with spatially-resolved observations. Most of the bursts 
in this study had rise times too short for observations of their progressive growth with 
the YLA, given its 10-s time resolution. Comparisons of the source areas from the VLA 
with burst rise times are possible and would provide an important test of these results. 
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CONCLUSION 


• Max ’91 study of impulsive hard X-ray/microwave bursts should focus on identify- 
ing the causative agent suggested herein 

• Combined images in hard X rays and microwaves ( e.g . GRID/VLA) should be used 
to search for the following phenomena. 

* Shock acceleration 

* Thermal conduction fronts 

* Traveling magnetic reconnection instabilities 

• These possible causative agents should reveal themselves via their distinctive rela- 

tionships with coronal magnetic structures 

* Quasi-perpendicular orientation of field lines relative to a radially expanding 

disturbance 

* Conductive growth of hot sources longitudinally along field lines 

* Longitudinal source expansion at speed near Alfven velocity 
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